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Abstract 

We review recent studies of light scalar meson properties by means of uni- 
tarization techniques, obtained from dispersion theory, and applied to the 
Chiral Perturbation Theory expansion. In particular, light scalars do not 
follow the N c dependence of qq states although a subdominant qq compo- 
nent may be observed to arise for the /o(600) around 1 GeV, where another 
qq multiplet is believed to exist. Finally, we present our preliminary results 
on the quark mass dependence of the / (600) and p(770) resonances. 

1 Introduction 

Light hadron spectroscopy lies outside the applicability range of QCD per- 
turbative calculations. Still, in this low energy region one can use Chiral 
Perturbation Theory (ChPT) [1] to obtain a model independent description 
of the dynamics of pions, kaons and etas. These particles are the Goldstone 
Bosons (GB) associated to the QCD spontaneous breaking of Chiral Symme- 
try and ChPT is built as a low energy expansion that contains those fields 
in the terms of a Lagrangian that respect all QCD symmetries, including its 
symmetry breaking pattern. The small quark masses of the three lightest fla- 
vors can be treated systematically within the perturbative chiral expansion 
and thus ChPT becomes a series in momenta and meson masses, generically 
0(p 2 /A 2 ). At lowest order there are no free parameters apart from masses 
and f n , the pion decay constant, that sets the scale A = 4Trf n . The chiral 
expansion can be renormalized order by order by absorbing the loop diver- 
gences in higher order counterterms, known as low energy constants (LEC), 
whose values depend on the specific QCD dynamics. That is, other theories 
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with spontaneous chiral symmetry breaking at the same scale will have the 
same leading order, but will differ in the values of the LEC. 

The renormalized LEC have to be determined from experiment, since they 
cannot be calculated from perturbative QCD. However, thanks to the fact 
that ChPT has the same symmetries than QCD and that it should couple 
to different kind of currents in the same way, it is still possible to determine 
in a model independent way how the constants that appear in ChPT, and 
therefore the observables, depend on some QCD parameters. This is indeed 
the case of the leading dependence on the number of colors N c and of the 
dependence on the quark masses. 



2 Unitarization and dispersion theory 

The unitarity of the S matrix implies that, for physical values of s, partial 
waves t IJ of definite isospin I and angular momenta J for elastic meson- 
meson scattering should satisfy 

lmt IJ = a\t IJ \ 2 Im = -a (1) 

where a = 2p/y / i, and p is the CM momenta of the two mesons. Note that 
unitarity implies that \t IJ \ < 1/cr, and a strong interaction is characterized 
precisely by the saturation of this unitarity bound. 

However, partial waves are obtained within ChPT as a low energy expan- 
sion t ~ t 2 + 1 4 + 1 6 + . . . , (To simplify the notation, from now on we will drop 
the I. J indices.) where t 2 k = 0(p/(47r/ 7r )) 2fe , and thus they cannot satisfy 
unitarity exactly, but just perturbatively, i.e: 

lmt 2 = 0, Imt 4 = atl, etc... (2) 

Unitarization methods extend the ChPT series to high energies by using 
the fact, remarked in eq.([T]), that the imaginary part of the inverse amplitude 
is known exactly. Hence, we can approximate the real part of Ret -1 ~ 
t 2 2 {t 2 + Ret 4 + ...) with ChPT, and find that 



Ret -1 — ia 1 — t i /t 2 

This is known as the one-channel Inverse Amplitude Method [2, 3] . A usual 
complaint is that the ChPT series is only valid at low energies, and there is 
no reason to use it beyond that regime. 

However, not only the complete t, but also the inverse amplitude 1/t 
and the ChPT series to next to leading order (NLO) and beyond, have an 
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analytic structure with a "physical cut" extending from threshold to oo and 
a "left cut" from — oo up to the highest value of s allowed in the crossed 
channel (s = if the two mesons are identical). It is then possible to write 
the following dispersion relation [3] for t& 

, 2 s 3 f°° Imt 4 (s')ds' r „. , 
U = b + b lS + b 2 s 2 + - / 4V ; + LC(U). 4 

where "LC" stands for a similar integral over the left cut and we have three 
subtractions to ensure convergence. A similar dispersion relation can be writ- 
ten for the function G = t\/t, by simply replacing i 4 by G and changing the 
name of the subtraction constants. Since £2 is real, these two functions have 
opposite imaginary parts on the physical cut so that their contributions from 
the "physical cut" integral are exactly opposite. Their subtractions constants 
correspond to the value of these functions at s = where it is perfectly jus- 
tified to use the ChPT expansion. And finally, their contributions from the 
left cut are also opposite from each other, but this time only up to the NNLO 
ChPT terms. This guarantees a good description of the integrand only at 
low energies, but that is precisely the region that has been weighted by the 
three subtractions. 

Therefore, note that the IAM is actually exact on the physical region and 
only uses the ChPT approximations for the subtractions constants at s — 0, 
where the use of ChPT is totally justified, and for the left cut, where the use 
of ChPT might not be justified at large \s'\, but the influence of this region 
is dumped by the subtractions. The use of the IAM is even more justified 
if used sufficiently far from this left cut, since this has an additional (s — s') 
suppression. In summary, there are no model dependencies in the approach, 
but just approximations to a given order in ChPT, and indeed the IAM can 
be derived not only at NLO but also at NNLO. 

Remarkably, the simple formula of the IAM, eq.pl), is able to describe 
the 7T7T and ttK scattering data not only at low energies, where it reproduces 
the ChPT series, but also in the resonance region. This is done with values 
of the ChPT parameters that are compatible with the values obtained within 
standard ChPT. 

In addition, the IAM generates the poles [3,4] associated to the resonances 
in the second Riemann sheet. This is of relevance since, in particular, the 
scalar resonances are the subject of intense debate that has been lasting for 
several decades, and as we have seen, the IAM is able to generate their poles 
from first principles like unitarity, analyticity and the QCD chiral symmetry 
breaking, without introducing these resonances by hand. Thus, we can study, 
without any a priory assumption, the nature of these states as follows from 
first principles and QCD. 
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3 N c behavior 

The QCD 1/N C expansion [5] allows for a clear identification of a qq res- 
onance, since it becomes a bound state, whose width follows an 0(1/N C ) 
behavior, whereas its mass should behave as 0(1). For our purposes, the rel- 
evant observation is that the leading l/N c behavior of the ChPT constants is 
known in a model independent way. Thus, in order to know the leading N c 
behavior of the resonances generated with the IAM, we just have to change 
the ChPT parameters according to their established N c scaling properties. 
For instance, the pion decay constant scales as 0(\/AQ, and we thus substi- 
tute f n by f w a/A^ c /3. Similar replacements, but according to their respective 
N c scaling, have to be done with all ChPT parameters. 

This procedure [4, 6] was first applied to the coupled channel IAM [4, 7] 
using the one loop SU(3) ChPT amplitudes, that also include kaons and 
etas, and the result was that the light vector resonances p(770) and if* (982) 
followed remarkably well the expected behavior of qq states. In contrast, 
the members of the light scalar nonet, namely, the sigma or / (600), the 
kappa or K (800), as well as the /o(980) and ao(980) resonances, all showed 
a behavior at odds with that of qq states (Only the a (980) could display 
a qq behavior in a limited corner of parameter space). It follows that the 
dominant component of light scalar mesons does not have a qq nature. Other 
large N c analysis have reached similar conclusions [8]. 

Very recently [9] we have extended the analysis to two loops in SU(2) 
ChPT, using the IAM one-channel formalism just described above to NNLO. 
In addition, we have developed a quantitative measure of how close a reso- 
nance is to a qq behavior: Taking into account subleading orders, to consider 
a resonance a qq state, it is enough that 

Ml q = M ( 1 + — V = — (i + ^L], (5) 

Nc V nJ n < n c \ n c J' k j 

were M and T are unknown but N c independent and the subleading terms 
have been gathered in €m, er — 1, since we expect generically 30% uncer- 
tainties at N c = 3. Thus, for a qq state, the expected M^ c and rV c can be 
obtained from those at N c — 1 generated by the IAM, as follows 

1 1 



1 + e M 



„ r^-! (iv c - 1) 



N r N c -1 

, 1 I 

1 + e r 



= M Nc _, + AM" (6) 



N r N r ~l 



N r + ^ 



Note the qq index for quantities obtained assuming a qq behavior. The reason 
to write the iV c values from those at N c — 1 is to be able to calculate from 
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what N c a resonance starts behaving as a qq, which is of interest to search for 
subdominant qq components. Thus, we can define an averaged x| g to check 
how close a resonance is to a qq behavior, as follows: 



When this quantity is smaller than one, it indicates a qq behavior, whereas 
a larger value indicates that it does not behave predominantly as such. Note 
also that imposing the minimization of the x 2 m eq.([7]) we could try to force 
a given qq behavior for a given resonance when fitting data. 

When evaluating eq.([7j) above, one has to be careful not to consider too 
large N c values. The reason is that, after all, we are interested in the physical 
state at N c — 3 and, most likely, states are a mixture of different components 
with different N c behavior. By allowing for too large values of N c we could 
be altering too radically the nature of the state, and, since qq states are 
expected to survive in the large iV c limit, whereas other kind of states, like 
tetraquarks, glueballs, etc... do not, even insignificant admixtures of qq at 
N c = 3 could become dominant for sufficiently large N c . Thus, the most 
relevant information will come from N c not too far from N c = 3 and we 
consider N c values smaller than one order of magnitude its physical value of 
3, let us say N c < 20. 

In this respect, J. A. Oiler raised an interesting concern at the end of 
my talk in this conference, about the absence of r/(980) in our calculations. 
Certainly, the r( mass is due to the Ua(1) anomaly, and decreases like l/y/N~ c . 
If we were to consider very large N c , this particle should be definitely included 
in our calculations, because it would become the most relevant degree of 
freedom of QCD at very large N c . However, if we limit ourselves to, say 
N c = 20, its mass would be 950 MeV x a/3/20 ~ 370 MeV, so, it is still much 
heavier than the pions that can still be considered as the only low energy 
degrees of freedom of the theory (Furthermore, its contribution would be 
similar to that of kaon loops, although there are both neutral and charged 
kaon loops, whereas we only have one 77'). This is an additional reason why 
it would be just wrong to consider too large values of N c to draw conclusions 
about the nature of the sigma. 

Therefore, with the measure in eq. ((TJ) and n = 20, we are able to quantify 
the deviation of the /o(600) from the qq behavior: at NLO, even in the most 
favorable case when we try to impose that it behaves as a qq, the data fit 
yields a % q = 125 for the / (600). 

When using the NNLO (two loop ChPT calculation), we have many more 
ChPT parameters, which are not well known, and give a great deal of free- 
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dom. For the NNLO low energy constants (LECS) we thus use standard 
estimates [10] with a 100% uncertainty. Still, there are large correlations 
and some very weak dependences on some of these parameters that can be 
driven far from their standard values for negligible improvements in the x 2 
of the data fit. For that reason we stabilize the values of the LECS, impos- 
ing also the minimization of a Xlecs together with that of the data. Also, 
using eq.([7|), we impose the p(770), which is a well established qq state, to 
behave as such. Thus, at two loops, we find that the / (600) comes out with 
X^ q = 4 in the most favorable case when we try to impose it to behave as a 
qq. Even relaxing the p qq behavior, we still get x| Q = 3.5 for the / (600). In 
conclusion, the two loop I AM confirms once again that the / (600) does not 
behave predominantly as a qq state, whereas that behavior is nicely followed 
by the p, whose x% q < 0-35 at NLO and x§ 9 = 0.93 at NNLO. 

We show in Fig{T] the NNLO behavior of the p and /o(600) mass and 
width. In the top figure, we show the N c = 3 behavior of the p. The full 
dots represent the values, for different N c , of its "pole mass" , M, whereas the 
empty dots represent its "pole width", T, both normalized to their physical 
values at N c = 3. It can be clearly seen that, already for very low N c , M/M 3 
starts behaving as 0(1) and r/r 3 as 0(1/N C ), as expected for a qq state. 
As commented above, the behavior shown in this plot yields x^ q = 0.93 for 
the p. In contrast, in the bottom figure, and with the same conventions, we 
show the N c behavior of the a or / (600). This time we are also imposing 
in the fit minimization that it should behave as a qq state. Obviously it 
does not, at least until N c ~ 8, and its \% q — 4. Despite the a is still not 
behaving predominantly as a qq state, the price due to trying to impose such 
behavior is that the fit suffers a clear deterioration, since the data now has 
a x 2 /dof = 1.5 (compared to 1.1 before) and the p now has x| 9 = 1.3, thus 
with a much worse qq behavior. 

Remarkably, it is also clear that the a or /o(600) follows a qq behavior 
for iV c > 8 or 10. This suggests the existence of a subdominant qq behavior 
in the a or / (600) that originates at a mass around twice that of the sigma 
~ lGeV, where a qq nonet is usually located. This is in good agreement with 
the emerging view that there could be two scalr nonets, a non qq one below 
1 GeV a a qq above [11]. 

4 Quark mass dependence 

We now present our very preliminary work [13] on the quark mass dependence 
of the mass and widths of the a and p mesons, which could be of interest to 
compare with lattice studies, where the small physical masses of light quarks 
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Real and Imaginary part of p pole position 

Imposing p to be qq 
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Real and Imaginary part of o pole position 

Imposing p and o to be qq 
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Figure 1: Mass and width dependence on N c for the p (Top) and a (Bottom) 
resonances. We have normalized all M and T to their respective physical 
values at N c = 3. Note that, already for N c close to 3, the p mass and width 
follow nicely the respective 0(1) and 0(1/N C ) behavior expected for a qq 
state. In contrast, the a does not behave predominantly as such, but only 
develops, above N c = 8 or 10, a subdominant qq behavior with a mass ~ 1 
GeV. 
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are hard to implement. 

With the IAM we can also study the quark mass dependence of the light 
resonances by changing the meson masses in the amplitudes, which is equiv- 
alent to change quark masses, since, to leading order ~ m q + .... We 
have done this in SU(2) Unitarized ChPT at NLO. In tctt scattering at NLO 
there appear four LECs h, - ■ ■ ,U- Phase shifts are almost insensitive to 
and 11, for which we take the values in [1]: Z3 = 0.8 ± 3.8, l\ = 6.2 ± 5.7, 
evaluated at a renormalization scale of p = 770 MeV. In contrast, we make 
a data fit for l\ and l r 2 , finding l\ = -3.7 ± 0.2, l r 2 = 5.0 ± 0.4, in fairly 
good agreement with standard values. Finally, when changing pion masses 
we have to take into account that amplitudes are customarily written [1] in 
terms of the p independent LECs I [1] and the physical pion decay constant 

fw = fo ( 1 + iQ„2fi '4 + • • • ) that depend explicitly on the pion mass. 

We have taken two criteria to set the applicability limit of our method, 
that is, the maximum value of the pion mass we can use. First, we do not 
want to spoil the chiral expansion, and second, we do not want the two-pion 
threshold to reach the KK threshold. Taking into account that SU(3) ChPT 
works well with a kaon mass of ~495 MeV, and that, according to NLO 
ChPT, if we set M n ~ 500 MeV the kaon mass becomes ~ 600 MeV, this 
means that for 500 MeV pions, hit scattering is still elastic for about 200 MeV 
above threshold. Hence the above criteria impose an applicability bound of 
M,,. ~ 500 MeV. To go beyond that we would need a coupled channel SU(3) 
formalism. 

We show in FigfS] how the p and /o(600) poles in the second Riemman 
sheet move as M T changes. Both the p and a mass increase with the pion 
mass, but that of the a grows faster. In addition, both widths decrease, 
partly due to phase space reduction (the two-pion threshold grows faster 
than both resonance masses). When the two-pion threshold reaches the a 
mass, its pole remains for a short while on the second sheet with a non-zero 
width but quickly reaches the real axis where it meets its conjugate partner 
from the upper plane and splits again into two poles corresponding to virtual 
bound states located on the real axis below threshold. As the pion mass 
keeps increasing, one of those '"virtual state"' poles moves toward threshold 
and jumps into the first sheet, whereas the other one remains in the second 
sheet. Although, of course, this happens for very large masses, such an 
analytic structure, with two very asymmetric poles in different sheets of an 
angular momentum zero partial wave, may be a signal of molecular structure, 
as discussed by M. Pennington in this conference. 

Finally, as the pion mass increases, the p pole moves toward the real 
axis and just when threshold reaches its mass it jumps into the real axis on 
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Figure 2: p and er complex plane pole movement with increasing pion mass. 
To ease the comparison of the pole position relative to the two-pion threshold 
we normalize by the pion mass that is changing. Note how the sigma pole 
moves toward the real axis below threshold where it splits in two virtual 
states, whereas the p pole just moves toward threshold. 



the first sheet, thus becoming a traditional bound state, while its conjugate 
partner remains on the second sheet practically at the very same position as 
the one in the first. 

A publication with further details is in preparation [13] including results 
of the / (600) and p(770) mass and width evolution with the pion mass 
as well as a comparison with other works and lattice results. Estimates of 
uncertainties and possibly an extension to the SU(3) coupled channel case 
are presently in progress. 
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